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DERIVATION  OF  MOS-UNIQUE  CONTINUATION  RATES 


SUMMARY 

To  manage  Army  personnel  effectively,  accurate  predictions  of 
the  future  personnel  status  of  each  military  occupational  speci- 
alty (MOS)  are  essential.  The  Army  has  developed  several  determi- 
nistic models  capable  of  projecting  personnel  strengths  by  MOS  up 
to  4 years  In  time.  In  making  projections,  the  models  require 
continuation  rates  (probabilities  that  individuals  will  remain  in 
service)  for  each  year  of  service.  The  continuation  rates  used  by 
these  models  have  been  derived  from  the  total  Army  population; 
thus,  predictions  have  been  accurate  only  at  the  total  force  level 
and  not  at  the  Individual  MOS  level  as  required.  However,  accu- 
rate predictions  could  be  made  at  the  MOS  level  if  the  continua- 
tion rates  were  computed  at  the  MOS  level. 

Derivation  of  continuation  rates  at  the  MOS  level  of  detail  is 
complicated  by  two  factors:  (1)  promotions  involving  a change  in 
MOS,  and  (2)  MOS  reclassifications  within  grade.  The  latter  is 
normally  used  to  balance  MOS  in  which  assets  greatly  exceed  or 
fall  short  of  authorizations.  In  the  methodology  develop. d,  the 
effects  of  the  two  factors  are  accounted  for  by  using  a directed 
graph  to  represent  the  total  force,  where  the  nodes  represent 
populations  for  each  MOS/Grade  pair  and  the  arcs  quantitatively 
represent  personnel  flows  between  nodes. 

To  develop  pure  continuation  rates  at  the  MOS  level,  the  Back- 
ward Star  Traversal  Algorithm  is  used  to  remove  the  promotion  and 
reclassification  effects.  Convex  combinations  of  the  pure  rates 
are  then  computed  for  each  node  of  the  graph  according  to  the  com- 
position of  the  node.  The  convex  combinations  are  used  as  esti- 
mates of  the  continuation  rates  by  MOS  and  grade  for  each  year  of 
service. 
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DERIVATION  OF  MOS-UNigUE  CONTINUATION  RATFS 

CHAPTER  1 
INTRODUCTION 


1-1.  BACKGROUND.  Effective  personnel  management  requires  the 
Army  to  be  able  to  predict  accurately  personnel  shortages  by  Mili- 
tary Occupational  Specialty  (MOS).  There  are  over  400  enlisted 
MOS  in  the  Army  and  the  skills  related  to  each  MOS  tend  to  be  nar- 
rowly defined.  The  skill  definitions  in  many  cases  are  directly 
related  to  specific  grades  and,  as  a result,  promotions  frequently 
involve  changes  in  MOS.  An  added  complication  is  introduced  by 
MOS  reclassifications  within  grade  which  results  from  attempts  to 
balance  chronic  surpluses  and  shortages. 

1-2.  PROBLEM.  The  Army  maintains  data  on  the  historical  popula- 
tions of  MOS  by  both  grade  and  year  of  service;  however,  in  the 
past  the  Army  has  been  unable  to  derive  MOS-unique  continuation 
rates*  from  these  data  because  of  the  changes  of  MOS.  The  diffi- 
culty has  been  that  true  losses  to  an  MOS  in  the  Army  could  not  be 
distinguished  from  apparent  losses  caused  by  changes  of  MOS.  In 
the  absence  of  MOS-unique  continuation  rates,  average  rates  based 
on  the  total  Army  population  have  been  applied  at  the  MOS  level  to 
project  the  force  over  time.  This  has  caused  predictions  to  be 
accurate  only  at  an  aggregated  level  and  not  at  MOS  level  as  re- 
quired. The  Personnel  Inventory  Analysis/Year  of  Service  (PIA/ 
YOS)  model-  is  used  by  the  Army  to  project  the  force  over  time. 
The  model  requires  two  key  sets  of  input  data:  continuation  rates 

by  MOS/grade/year  of  service  and  "feeder  patterns"  which  specify 
the  direction  and  relative  magnitude  of  personnel  flows.  The  mo- 
del first  computes  updated  populations  by  grade  and  MOS  using  the 
continuation  rates  and  then  simulates  promotion  and  reclassifica- 
tion according  to  the  feeder  patterns.  Since  promotion  is  modeled 
explicitly,  the  continuation  rates  used  must  reflect  only  true 
losses  to  the  Army.  The  average  continuation  rates  based  on  the 
total  Army  population  as  previously  used  do  not  accurately  reflect 
true  losses  at  the  MOS  level. 

% 

— 

I 

‘Continuation  rates  are  defined  as  probabilities:  the  contin- 

uation rate  forbear  of  service  1 is  the  probability  that  an  indi- 
vidual in  his  v year  of  service  will  remain  in  service  at  least 
one  more  year. 
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1-3.  PURPOSE.  This  paper  describes  a technique  by  which  more 
precise  MOS-unique  continuation  rates  needed  by  the  PIA/YOS  model 
can  be  derived.  The  PIA/YOS  feeder  patterns,  which  also  serve  as 
Input  to  the  PIA/YOS  model,  are  first  restructured  as  a directed 
graph.  The  graph  is  then  used  in  conjunction  with  historical  data 
on  MOS  populations  to  derive  continuation  rates  which  are  unique 
by  MOS,  grade  and  year  of  service  and  yet  exclude  the  effects  of 
promotion.  These  rates  are  then  substituted  for  the  previously 
defined  rates  based  on  the  total  Army  populations  and  used  as  In- 
put to  the  PIA/YOS  model. 
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CHAPTER  2 

DERIVATION  OF  MOS-UNIQUE  CONTINUATION  RATES 


2.1  GENERAL.  In  the  directed  graph  formulation,  each  combination 
of  an  MOS  and  grade  Is  considered  as  a node  and  both  promotion  and 
reclassification  are  represented  as  arcs.  Before  generalizing 
this  concept,  consider  the  system  shown  in  Figure  2-1  which  has 
only  two  MOS,  A and  B,  and  two  grades,  El  and  E2. 


Figure  2-1.  A Notional  Example  Illustrating  Node  and  Arc 

Information 
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In  the  directed  graph,  the  arcs  represent  potential  flows  between 
nodes;  however,  the  amount  of  the  flow  is  determined  by  the  un- 
filled requirements  at  the  receiving  node.  In  Figure  2-1,  some 
unfilled  requirement  may  exist  at  node  N2  after  attrition.  If  a 
requirement  does  exist,  80%  of  that  requirement  will  be  filled 
from  node  N1  and  20%  from  node  N4.  Calculation  of  continuation 
rates  for  even  the  simplified  example  Is  not  straightforward.  The 
sub-populations  at  nodes  N1  and  N3  can  be  categorized  easily  since 
these  are  entry  level  grades.  Nodes  N2  and  N4  are  not  so  easily 
described.  Although  the  sub-population  of  N2  currently  holds  MOS 
A and  grade  E2,  at  least  part  of  that  sub-population  consists  of 
people  who  originally  held  MOS  B.  To  further  complicate  matters, 
some  of  the  people  serving  at  node  N2  who  originated  in  MOS  A have 
also  served  In  MOS  B having  moved  on  the  arcs  from  N1  to  N4  and 
then  N4  to  N2.  The  "Backward  Star  Traversal  Algorithm"  1BSTAR) 
developed  In  this  paper  uses  the  negative  star  operator  — to  com- 
pensate for  the  effects  of  MOS  flows  and  produce  the  required  MOS- 
unique  continuation  rates. 

2-2.  DESCRIBING  THE  GRAPH.  To  generalize  the  example  of  Figure 
2-1,  a method  for  representing  total  Army  enlisted  personnel  flows 
Is  described  using  a graph  G.  Sub-graphs  of  G corresponding  to 
each  Career  Management  Field  (CMF)  are  defined  and  a technique  Is 
presented  for  traversing  the  sub-graphs  to  compute  "Distribution 
Vector  Sets"  for  the  nodes  of  the  sub-graphs.  In  this  paper,  n 
represents  the  number  of  MOS  and  £ will  signify  the  number  of 
grades.  Where  m.  denotes  MOS  1 and  g,  denotes  grade  j,  the  graph 
G with  a node  set  V and  an  arc  set  A Is  defined  as  follows: 


(1)  G * (V,A),  where 

(2)  V * (vk  = (m1 ,9j ) | 1 = l,2,...n;  j = 1,2,...,*} 

for  k * 1,2,...  I VI  where  I VI  Is  the  cardinality 
of  the  set  V 

(3)  A = (ax  * (vs,vt,pst)  | s Mi  s - 1,2 |V|  ; t = 1,2 |V| 

0 - p$t  51}  for  x « 1,2,...,  |A| 


t: — t - 


_ 
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Because  grades  El  through  E3  are  frequently  considered  together, 
the  full  range  of  enlisted  grades  (E1-E9)  are  modelled  with  i 
equal  to  7.  The  number  of  MOS  could  vary  from  a low  of  400  to  an 
upper  bound  of  600  depending  on  the  time  frame  of  the  data.  At- 
tempting to  solve  the  problem  at  the  total  Army  level  presented 
two  major  computer-related  problems.  First,  computational  time  is 
directly  proportional  to  the  number  of  MOS  and  it  was  estimated 
that  a solution  of  a complete  graph  would  require  from  4 to  6 CPU 
hours.  Secondly,  core  requirements  would  exceed  the  limits  of 
most  computing  facilities  and  could  require  segmentation  of  the 
problem.  Fortunately,  these  problems  can  be  avoided  because  MOS 
are  currently  managed  and  controlled  within  disjoint  Career  Man- 
agement Fields.  If  X is  the  number  of  CMF's,  sub-graphs  of  G 
can  be  defined  in  the  following  manner: 

(4)  ^ - (V*,) 

X 

(5)  U = G;Gj  n Gj  = <p  v 1 ,j  3 1 t j 
i = l 

where  <l>  is  the  null  set 

X 

(6)  U = V;  ^nVj  = *vi,j  s i t j 

i=l 

X 

(7)  U ffj  = A;  frjOA.  = 4>  \/  i , j 3 1 f j 
1=1  J 


Neither  G nor  the  sub-graphs  G-j  are  directed  graphs  as  defined. 
However,  In  order  to  conform  to  the  feeder  patterns,  G and  there- 
fore Gi  must  be  considered  directed  graphs  of  a special  nature. 

If  v$  * (m.i , g$)  and  v^  = (m-j,  gt),  then  the  feeder  patterns  re- 
quire that: 

(8)  gt!gs  V a *(Wpst)  e A 

(9)  Circuits  are  not  allowed 


Constraint  (8)  requires  a left-to-right  and  up-or-down  ordering 
corresponding  to  promotion  or  within  grade  MOS  reclassification 
and  (9)  Is  needed  to  Insure  that  the  algorithm  will  terminate. 
Figure  2-2  Illustrates  the  basic  arc  types  allowed  by  the  con- 
straints. 
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vx  = (mk’9j+i) 


Figure  2-2.  Arc  Types  in  the  Directed  Graph 

In  Figure  2-2,  arc  a?  is  a promotion  arc  since  the  MOS  remains  the 
same,  but  the  grade  increases.  Arc  83  is  a within  grade  reclassi- 
fication arc  since  the  grade  remains  the  same  and  the  MOS  changes. 
Arc  ai  Is  a promotion  with  reclassification  arc  since  both  the 
grade  and  MOS  change.  To  traverse  the  sub-graphs,  the  BACKWARD 
STAR  OPERATOR  Is  defined:  

(10)  3{vt)  = (ae^  | a = (vx,vt,pxt)  V vxeVj} 

Thus  3{v^}  Identifies  all  arcs  which  terminate  at  node  v*.  In 
Figure  2-2,  3{ vx>  = {a^.aj}  and  3{vt)  r (82) 

2-3.  DESCRIBING  PRIOR  MOS  FLOWS.  The  set  D =< Di , D2,  . . . , Dk } 
where  k = |Vj|  is  defined  to  be  a "Distribution  Vector  Set,"  where 
the  elements  in  D are  In  one-to-one  correspondence  with  the  elements 
(nodes)  In  D Is  an  n-tuple  of  the  form  Dx  = (d],  do,  ....  dn) 
where  n Is  now  the  number  of  MOS  In  the  CMF.  The  purpose  of  the 
distribution  vector  is  to  summarize  Information  on  prior  flows  to  a 
given  node.  The  entries  In  Dx  should  contain  the  fraction  of  node 
vx  that  originated  in  MOS  1,  flOS  2,  ...,  MOS  n.  Thus  o^d|<5l 
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for  k = 1,2,  ....  h and  2 d.  = 1 . For  computational  ease,  an  n 

k=1  K 

by  £ matrix  r is  established  and  the  node  set  V is  mapped  into  T 
such  that 

1.  if  v = (mi,g.j)eV 


0,  otherwise 


The  n by  i by  n tensor  A is  also  established  and  D is  mapped 
into  A.  This  mapping  process  is  shown  in  Figure  2-3. 
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Since  the  elements  In  D are  In  one-to-one  correspondence  with  the 
nodes  In  V,  It  follows  that  an  automorphism  exists  which  maps  A 
Into  T with  respect  to  the  first  two  Indlcies.  The  tensor  A will 
be  zero  filled  Initially,  and  I'  will  be  used  to  determine  the  ap- 
propriate entries  for  A.  since  the  flows  between  nodes  are  left- 
to-right  or  up-or-down,  T Is  processed  by  column  from  left-to- 
rlght,  that  is,  column  j of  r Is  completely  processed  before  pro- 
cessing column  j+1.  When  processing  Is  completed,  A<M,k  contains 
the  percent  of  MOS  1,  grade  j that  originated  in  MOS'k')  Since  the 
entries  In  fare  surrogates  for  the  node  set,  an  entry  Tj  i 
cannot  be  processed  until  every  node  which  sends  flow  to  node  v = 
(m^  g^  has  been  processed.  To  record  such  processing,  the  STATE 
OPERATOR  S on  f Is  defined  such  that: 


s(  ri.j) 


1*  If  r ] ' j has  been  completely  processed 
1,  If  3{ v=(m1  ,gj ) } = <J> 

0,  otherwise 


To  process  the  jib  column  of  r , the  set 


Tj  * {1  I r itj  * i;S(r1tj)  ■ 0} 

Is  defined.  Thus,  T.  Is  the  set  of  entries  In  the  j—  column  of  T 
. Next,  backward  star  sets  for  each  element  of  T,  are  defined  as 
follows:  J 


(11)  B1,j  = <Kvx  “ (mi ,g j ) > y vxeV3leTj 


Using  the  backward  star  sets  formed  In  (11),  the  Individual  ele- 
ments In  Tj  can  be  processed  If,  and  only  If,  the  state  operator  S 
for  all  source  nodes  In  the  backward  star  set  are  equal  to  1;  that 
Is  the  1—  element  In  Tj  can  be  processed  Iff 

O2)  s(ra,b)  * 1 V vx  « (ma»9b)eBl  J 

If  (12)  Is  satisfied  for  the  1th  element  In  T«,  then  A<  < l Is 
calculated  for  each  k as  follows:  J 

*13)  A1.j.k  * A1,j,k  + Pst  Ax,y,k 

Vx,y3Vs  * (n»Xt9y)  e B1,J  for  k " 1 .2 n 
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When  (13)  has  been  completed  for  all  source  nodes  and  for  all  val- 
ues of  k.r^j  has  been  processed.  The  state  operator  S(  r i , i ) Is  set 
to  one  and  the  set  Ti  Is  updated:  T,  = T.  - ij  . If  (12)  Is  not 
satisfied,  that  element  of  T,  ii  bypdssedJand  the  next  element  of 
the  Tj  set  Is  processed.  IfJafter  sequentially  examining  and  pro- 
cessing the  elements  In  T.,  the  set  Is  null,  then  column  j has 
been  completely  processed':  If  T,  Is  not  null,  the  algorithm 

starts  again  with  the  first  element  In  T.  and  continues  processing 
until  T,  Is  null.  Note  that  because  of  constraints  (8)  and  (9), 
every  pass  through  T,  reduces  the  cardinality  of  the  set  by  at 
least  1.  Since  IT .1  J 1 s finite,  T,  must  become  null  after  a finite 
number  of  passes.  J J 

2-4.  Computing  Interim  MOS  populations.  The  a tensor  has  estab- 
lished an  audit  trail  of  movement  of  personnel  between  MOS  and  is 
now  used  to  form  Interim  MOS  as  If  no  change. of  MOS  had  occurred. 

For  every  MOS,  define  a tensor  P such  that  P1  Is  the  strength 
of  MOS  1 for  year  of  service  y,  grade  g and  d(l8fl8ar  quarter  . A 
tensor  I Is  also  defined  dimensioned  by  year  of  service,  MOS  and 
quarter  such  that  I Is  the  strength  of  interim  MOS  m In  year 
of  service  y In  quart£fMq.  The  entries  In  I are  computed  as  fol- 
lows: 


(14)  I 


y.m.q 


n 

= 2 
1=1 


1 

2 

j=l 


t‘WW] 


Note  that  (14)  reallocates  to  Interim  MOS  m everyone  that  origi- 
nated In  MOS  m at  the  entry  level.  Having  compensated  for  move- 
ment between  MOS,  Interim  MOS  continuation  rates  are  computed. 

2-5.  COMPUTING  INTERIM  MOS  CONTINUATION  RATES.  Equations  (15) 
and  (16)  Illustrate  two  different  methods  for  computing  these 
rates.  Let  T m,y  be  the  continuation  rate  for  Interim  MOS  m In 
the  yth  year  of  service.  One  computational  method  Is: 

8 

q*5 

05)  *ln,y  = 
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The  other  method  is: 


<16>  Vy 


4 T 

1 -r  1y+1,m,q+4 

4 L 

qa1  ^.m.q 


Both  methods  are  heuristic;  (lb)  Is  a weighted  moving  average  and 
(16)  Isa  simple  moving  average.  Both  use  eight  quarters  of  data 
and  both  are  candidates  for  use;  however,  additional  testing  Is 
required  before  a final  selection  Is  made. 

2-6.  MOS-UNIQUE  CONTINUATION  RATES.  With  the  Interim  continua- 
tion rates,  Acan  be  used  to  form  convex  combinations  of  the  In- 
terim MOS  rates  yielding  continuation  rates  which  are  MOS-unlque 
at  the  grade  and  year  of  service  level  of  detail.  The  computa- 
tional technique  Is  given  by  (17)  where  Ri  . Is  the  rate  for  MOS 
k 1,  grade  j and  year  of  service  y. 

(.A1.J.kMk,y)»  iff  r1,j  = 1 

0,  otherwise 


(17)  R 
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A measure  of  the  validity  of  the  methodology  can  be  obtained  by 
applying  the  rates  to  the  data  on  which  the  rates  are  based.  June 
1975  and  June  1976  data  were  used  to  compute  Army  wide  continuation 
rates  and  MOS-unlque  continuation  rates  for  the  maneuver  combat 


arms.  Both  rates  were  used  to  "predict"  the  1976  first  year  of 
service  population  from  the  June  1975  data  for  MOS  11 B,  light 
weapons  Infantryman.  The  prediction  obtained  using  MOS-unlque 
rates  was  within  0.27%  of  actual  experience  while  Army  wide  rates 
were  In  error  by  69.5%.  While  this  is  a biased  test.  It  appears 
that  the  use  of  MOS-unlque  rates  offers  a significant  advantage 
over  the  current  method  of  applying  Army  wide  rates  at  the  MOS 
level. 
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